Indication for the disappearance of reactor v e in the Double Chooz experiment 
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The Double Chooz Experiment presents an indication of reactor electron antineutrino dis- 
appearance consistent with neutrino oscillations. An observed-to-predicted ratio of events of 
0.944 ± 0.016 (stat) ± 0.040 (syst) was obtained in 101 days of running at the Chooz Nuclear Power 
Plant in France, with two 4.25 GWu reactors. The results were obtained from a single 10 m 3 fiducial 
volume detector located 1050 m from the two reactor cores. The reactor antineutrino flux prediction 
used the Bugey4 flux measurement after correction for differences in core composition. The deficit 
can be interpreted as an indication of a non-zero value of the still unmeasured neutrino mixing 
parameter sin 2 2#i3. Analyzing both the rate of the prompt positrons and their energy spectrum we 
find sin 2 26»i 3 = 0.086 ± 0.041 (stat) ±0.030 (syst), or, at 90% CL, 0.017 < sin 2 26»i 3 < 0.16. 

PACS numbers: 14.60. Pq,13.15.+g,25.30.Pt,95.55.Vj,28.41.Ak 
Keywords: neutrino oscillations, neutrino mixing, reactor 



We report first results of a search for a non-zero neu- 
trino oscillation |Ti mixing angle, #13, based on reactor 
antineutrino disappearance. This is the last of the three 
neutrino oscillation mixing angles @,Q for which only up- 
per limits 4, 5] are available. The size of #13 sets the re- 
quired sensitivity of long-baseline oscillation experiments 
attempting to measure CP violation in the neutrino sec- 
tor or the mass hierarchy. 

In reactor experiments [1, 0] addressing the disap- 
pearance of v e , 0i3 determines the survival probability 
of electron antineutrinos at the "atmospheric" squared- 
mass difference, Am 2 tm . This probability is given by: 

P surv « 1 - sin 2 26 13 sm 2 (l.267Am 2 atm L/E) , (1) 

where L is the distance from reactor to detector in meters 
and E the energy of the antineutrino in MeV. The full for- 
mula can be found in Ref. [J). Eq.[T] provides a direct way 
to measure #13 since the only additional input is the well 
measured value of | Am 2 t J = (2.32±H|}x 10" 3 eV 2 [|. 
Other running reactor experiments [9l llOj are using the 
same technique. 

Electron antineutrinos of < 9 MeV are produced by 
reactors and detected through inverse beta decay (IBD): 
9 e + p — > e + + n. Detectors based on hydrocarbon liquid 
scintillators provide the free proton targets. The IBD 
signature is a coincidence of a prompt positron signal 
followed by a delayed neutron capture. The v e energy, 
E Pc , is reconstructable from E promp t , the positron visible 
energy (E Pc S E prompt + 0.78 MeV). 

Recently, indications of non-zero #13 have been re- 
ported by two accelerator appearance experiments: 
T2K [llf and MINOS Q3. Global fits (see 
e.g. [la, [l4|) indicate central values in the range 
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0.05 < sin 2 2fli 3 < 0.10, accessible to the Double Chooz 
experiment [H, [l6[ . 

We present here our first results with a detector lo- 
cated ~ 1050 m from the two 4.25 GWth thermal power 
reactors of the Chooz Nuclear Power Plant and under a 
300 MWE rock overburden. The analysis is based on 101 
days of data including 16 days with one reactor off and 
one day with both reactors off. 

The antineutrino flux of each reactor depends on its 
thermal power and, for the four main fissioning isotopes, 
235 U, 239 Pu, 238 U, 241 Pu, their fraction of the total fuel 
content, their energy released per fission, and their fis- 
sion and capture cross-sections. The fission rates and as- 
sociated errors were evaluated using two predictive and 
complementary reactor simulation codes: MURE [l~7l [l~8l] 
and DRAGON 19]. This allowed a study of the sensi- 
tivity to the important reactor parameters (e.g.. ther- 
mal power, boron concentration, temperatures and den- 
sities). The quality of these simulations was evaluated 
through benchmarks (20| . and comparisons with Elec- 
tricite de France (EDF) assembly simulations. The max- 
imum discrepancies observed were included in the fission 
rate systematic error. 

MURE was used to develop a 3D simulation of the re- 
actor cores. EDF provided the information required to 
simulate the fission rates including initial burnups of as- 
semblies. To determine the inventories of each assembly 
composing the core at the startup of the data-taking cy- 
cle, assembly simulations were performed and the inven- 
tories at the given burnup computed. The energies per 
fission computed by Kopeikin [21[ and nuclear data eval- 
uated from the JEFF3.1 database [22j were used. The 
evolutions of the core simulations with time were per- 
formed using the thermal power and the boron concen- 
tration from the EDF database averaged over 48 h time 
steps, yielding the relative contributions to fissions of the 
four main isotopes. 

The associated antineutrino flux was computed using 
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the improved spectra from [23l, c onverted from the ILL 
reference electron spectra [2414261 ] . and the updated ab 
initio calculation of the 238 U spectrum [27j. The ILL 
spectra were measured after irradiating U or Pu for ~ 1 
day. Contributions from /3-decays with lifetimes longer 
than these irradiation times were accounted for as pre- 
scribed in (27j . 

The Double Chooz detector system (Figure [T]) consists 
of a main detector, an outer veto, and calibration devices. 
The main detector comprises four concentric cylindrical 
tanks filled with liquid scintillators or mineral oil. The 
innermost 8 mm thick transparent (UV to visible) acrylic 
vessel houses the 10 m 3 ^-target liquid, a mixture of n- 
dodecane, PXE, PPO, bis-MSB and 1 g gadolinium/1 as 
a beta-diketonate complex. The scintillator choice em- 
phasizes radiopurity and long term stability [28|. The 
i/- target volume is surrounded by the 7-catcher, a 55 cm 
thick Gd-free liquid scintillator layer in a second 12 mm 
thick acrylic vessel, used to detect 7-rays escaping from 
the {/-target. The light yield of the 7-catcher was chosen 
to provide identical photoelectron (pe) yield across these 
two layers [Hj]. Outside the 7-catcher is the buffer, a 
105 cm thick mineral oil layer. It shields from radioac- 
tivity of photomultipliers (PMTs) and of the surround- 
ing rock, and is one of the major improvements over the 
CHOOZ experiment 0. 390 10-inch PMTs [liHH are 
installed on the stainless steel buffer tank inner wall to 
collect light from the inner volumes. These three volumes 
and the PMTs constitute the inner detector (ID). 

Outside the ID, and optically separated from it, is a 
50 cm thick "inner veto" liquid scintillator (IV). It is 
equipped with 78 8-inch PMTs and functions as a cos- 



mic muon veto and as a shield to spallation neutrons pro- 
duced outside the detector. The detector is surrounded 
by 15 cm of demagnetized steel to suppress external 7- 
rays. The main detector is covered by an outer veto 
system (not used in this analysis). 

The readout is triggered by custom energy sum elec- 
tronics [HI- 35]. The ID PMTs are separated into two 
groups of 195 PMTs uniformly distributed throughout 
the volume and the PMT signals in each group are 
summed. The signals of the IV PMTs are also summed. 
If any of the three sums is above a set energy threshold, 
the detector is read out with 500 MHz flash- ADC elec- 
tronics [36H37I ] with customized firmware and a deadtime- 
free acquisition system. Upon each trigger, a 256 ns 
interval of the waveforms of both ID and IV signals is 
recorded. The low trigger rate (120 Hz) allowed the ID 
readout threshold to be set at 350 keV, well below the 
1.02 MeV minimum energy of an IBD positron, greatly 
reducing the threshold systematics. 

The experiment is calibrated by several methods. A 
multi- wave length LED-fiber light injection system (LI) 
produces fast light pulses illuminating the PMTs from 
fixed positions. Radio-isotopes 137 Cs, 68 Ge, 60 Co, and 
252 Cf were deployed in the target along the vertical sym- 
metry axis and, in the gamma catcher, through a rigid 
loop traversing the interior and passing along boundaries 
with the target and the buffer. The detector was mon- 
itored using spallation neutron captures on H and Gd, 
residual natural radioactivity, and daily LI runs. The 
stability of the peak energy of neutron captures on Gd 
in IBD candidates is shown in Figure [2j The energy re- 
sponse was found to be stable within 1% over time. 
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FIG. 1. A cross-sectional view of the Double Chooz detector 
system. 
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FIG. 2. The peak of the energy of neutron captures on Gd 
in IBD events(right scale) and its deviation from its average 
value(left scale) as a function of elapsed(calendar) day. 

The signature of IBD events is a delayed coincidence 
between a prompt positron energy deposition, E promp t, 
and a delayed energy deposition, Edoiay, due to the neu- 
tron capture on H or Gd within At e + n . The fiducial 
volume is constrained to the target vessel without posi- 
tion cuts by requiring a v e event to have a capture on 
Gd, identified by its emission of ~ 8 MeV in 7 rays. The 
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analysis compares the number and energy distribution of 
detected events to a prediction based on the reactor data. 

Energy measurements are based on the total charge, 
Qtoi, collected by the PMTs and corrected for gain vari- 
ations. The energy is reconstructed scaling Q tot by a con- 
stant, adjusted so that the energy of the gamma emitted 
following neutron capture on H reconstructs to 2.22 McV 
at the target center. This corresponds to ~ 200 pe/MeV. 
Our Monte Carlo (MC), based on GEANT4 0, is used 
to model the detector response and to calculate its accep- 
tance. It uses parameters for quenching [39(, absorption, 
re-emission, refraction, etc. determined from laboratory 
measurements of the detector liquids. Comparisons be- 
tween actual and simulated calibration data were used 
to develop a parametric function to correct the simula- 
tion, and to assess the uncertainties in the energy re- 
construction. The function is a product of two factors. 
One, dependent on energy, ranges from 0.97 to 1.05 for 
0.7-10.0 MeV. The other, dependent on position, ranges 
from 0.94 to 1.00 over the target volume. 

The following criteria are applied to select 
v e candidates. Triggers within a 1000 fis window 
following a cosmic muon crossing the IV or the ID 
(46 s _1 ) are rejected to limit spallation neutron and 
cosmogenic backgrounds. This requirement is followed 
by five selections: 1) a cut rejecting events caused by 
some sporadically glowing PMT bases, resulting in 
light localized to a few PMTs and spread out in time: 
Qmax/Qtot <0.09 (0.06) for the prompt (delayed) energy 
and rms(t s tart)< 40 ns, where Q m ax is the maximum 
charge recorded by a single PMT and rms(t s tart) is 
the standard deviation of the times of the first pulse 
on each PMT; 2) 0.7 MeV<E prompt <12.2 MeV; 3) 
6.0 MeV <E dc i ay <12.0 MeV; 4) 2 ^s < At e+n < 100 /is, 
where the lower cut eliminates correlated noise and 
the upper cut is determined by the ~ 30 /is capture 
time on Gd; 5) a multiplicity cut to reject correlated 
backgrounds defined as no additional valid trigger from 
100 /is preceding the prompt candidate to 400 /is after 
it. Applying selections (1-5) yields 4121 candidates or 
42.6 ± 0.7 events/day, uniformly distributed within the 
target, for an analysis live time of 96.8 days. 

Contributions from background events surviving these 
cuts have been estimated as follows. Uncorrelated coin- 
cidences result mainly from the random association of a 
prompt energy deposition due to radioactivity (7.6 s _1 ) 
and a later candidate neutron capture (~20/hour). This 
background is measured by applying selection cuts (1-5) 
but modifying selection (4) such that the 2 — 100 /is time 
window is shifted by 1000 fis relative to the prompt trig- 
ger. To improve the precision of this background mea- 
surement, 198 such windows, each shifted from the previ- 
ous one by 500 fj,s, were used, leading to 0.33±0.03 events 
per day. 

Fast neutrons induced by muons traversing the rock 
can interact in the target producing a recoil proton and, 
later, be captured, simulating an IBD event. We es- 
timate this rate to be 0.83 ± 0.38 events per day (in- 



cluding a contribution from stopping muons) by apply- 
ing cuts (1-5), but modifying selection (2) such that 
12.2 MeV <E prompt < 30 MeV, and then extrapolating to 
the signal region, assuming a flat energy spectrum. We 
account for an uncertainty in this extrapolation, and for 
the contribution of stopping muons, by including a shape 
error ranging up to ±70% of the flat extrapolation at 
lower energies. 

9 Li /3-n emitters are produced preferentially by en- 
ergetic muons. They were studied by searching for 
a triple delayed coincidence between a muon deposit- 
ing > 600 MeV in the detector and a D e -\ike pair of 
events, where the delay between the muon and prompt 
event is dictated by the 178 ms 9 Li halflife, which pre- 
cludes vetoing on all muons. Fitting the resulting time 
distribution with a flat component and an exponential 
with the 9 Li lifetime results in an estimated rate of 
2.3 ± 1.2 events/day. This rate is assigned the energy 
spectrum of the 9 Li decay branches. A shape uncertainty 
of up to 20% is introduced to account for uncertainties 
in some decay branches. 8 He is not considered since it 
is less abundantly produced [40]. The total background 
rate, 3.46 ± 1.26 d _1 , is summarized in TableQ] 

The overall background envelope is independently ver- 
ified by analyzing 22.5 hours of both-reactors-off data 
(< 0.3 residual v e events). Two D e candidates, with 
prompt energies of 4.8 MeV and 9.8 MeV, pass cuts (1- 
5). They were associated within 30 cm and 220 ms with 
the closest energetic muon, and are thus likely to be as- 
sociated with 9 Li. 



TABLE I. The breakdown of the estimated background rate. 
Additional shape uncertainties are described in the text. 



Background Rate/day Syst. 


Uncertainty (% of signal) 


Accidental 0.33 ± 0.03 


< 0.1 


Fast neutron 0.83 ± 0.38 


0.9 


9 Li 2.3 ±1.2 


2.8 



The following detector-related corrections and efficien- 
cies as well as their uncertainties were evaluated using 
the MC. The energy response introduces a 1.7% system- 
atic uncertainty determined from fits to calibration data. 
The number of free protons in the target scintillator, 
6.747 x 10 29 based on its weight measurement, has an 
uncertainty of 0.3%, originating from the knowledge of 
the scintillator hydrogen ratio. A dedicated simulation 



including molecular bond effects [41| indicates that the 
number of IBD events occurring in the gamma catcher 
with the neutron captured in the target (spill in) exceeds 
the number of events in the target with the neutron es- 
caping to the gamma catcher (spill out) by 1.4% ± 0.4%, 
0.8% lower than our standard MC prediction which was 
therefore reduced accordingly. Above the 700 keV anal- 
ysis threshold, the trigger efficiency is 100.0l~Q 4 %, as- 
sessed with a low threshold prescaled trigger. Calibra- 
tion data taken with the 252 Cf source were used to check 
the MC for any biases in the neutron selection criteria 
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and estimate their contributions to the systematic un- 
certainty. The fraction of neutron captures on Gd is 
found to be (86.0 ± 0.5)% near the center of the target, 
2.0% lower than the simulation prediction which was re- 
duced accordingly with a relative systematic uncertainty 
of 0.6%. The simulation reproduces the 96.5% efficiency 
of the Ai e + ra cut with an uncertainty of 0.5% and the 
94.5% fraction of neutron captures on Gd accepted by 
the 6.0 MeV cut with an uncertainty of 0.6%. The MC 
normalization was adjusted for the muon veto (—4.5%) 
and the multiplicity veto (—0.5%) dead-times. 



TABLE II. Contributions of the detector and reactor errors 
to the absolute normalization systematic uncertainty. 



Detector 




Reactor 




Energy response 


1.7% 


Bugey4 measurement 1.4% 


Edday Containment 


0.6% 


Fuel Composition 


0.9% 


Gd Fraction 


0.6% 


Thermal Power 


0.5% 




0.5% 


Reference Spectra 


0.5% 


Spill in/out 


0.4% 


Energy per Fission 


0.2% 


Trigger Efficiency 


0.4% 


IBD Cross Section 


0.2% 


Target H 


0.3% 


Baseline 


0.2% 


Total 


2.1 % 


Total 


1.8% 



The full covariance matrix of the emitted v e spectra 
was computed as prescribed in [27| . MURE pro- 
vided the fractions of fissions per isotope 235 U=48.8%, 
239 Pu=35.9%, 241 Pu=6.7%, and 238 U=8.7% and the fis- 
sion rate covariance matrix. The resulting relative uncer- 
tainties on the above fission fractions are ±3.3%, ±4%, 
±11.0% and ±6.5%, respectively. The error associated 
with the thermal power is ±0.46% at full power (42l. |43|. 
fully correlated between the two cores. 

To avoid being affected by possible very short base- 
line v e oscillations 0,13 |45||, we adopt the reactor 
v e spectrum of [H, |27|, but the global normalization 
is fixed by the Bugey4 rate measurement (4(| with its 
associated 1.4% uncertainty. A relative correction of 
(0.9 ±1.3%) of the Bugey4 value accounts for the dif- 
ference in core inventories. The IBD differential cross 
section is taken from [47j], using 881.5 ± 1.5 s [l| as the 
neutron lifetime. The systematic uncertainties are sum- 
marized in Table [TXJ The expected no-oscillation number 
of v e candidates is 4344 ± 165, including background. 

The measured daily rate of IBD candidates as a func- 
tion of the no-oscillation expected rate for different reac- 
tor power conditions is shown in Figure [3] The extrapola- 
tion to zero reactor power of the fit to the data (including 
the both-reactors-off) yields 3.2 ± 1.3 events per day, in 
excellent agreement with our background estimate and 
the both-reactors-off data. 

Our measurement can be expressed as an ob- 
served IBD cross section per fission, a^ c , a quan- 
tity which depends on the number of events ob- 
served, the number of target protons, the detector ef- 
ficiency, the number of fissions occurring during our 
measurement and the distance to the reactors, yielding 
af c = (5.383 ± 0.210) 10" 43 cm 2 /fission. The Bugey4 
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FIG. 3. Daily number of v R candidates as a function of the 
expected number of . The dashed line is a fit to the data, 
the band is the 90% C.L. of this fit. The dotted line is the 
expectation in the no-oscillation scenario. The triangle indi- 
cates the measurement with both reactors off. 



measurement, corrected to match our fractions of iso- 
topes quoted above, yields a cross section per fis- 
sion of (5.703 ± 0.108) 10~ 43 cm 2 /fission. The ratio of 
these two measurements is independent of any possi- 
ble very short baseline oscilations. (Without Bugey4 
normalization, the prediction, for our running con- 
ditions and using the reference spectra (23l . l27j . is 
(6.209 ± 0.170) 10" 43 cm 2 /fission). 

The ratio of observed to expected events is 
R DC = 0.944 ± 0.016 (stat) ± 0.040 (syst), correspond- 
ing to sin 2 26>i 3 = 0.104 ± 0.030 (stat) ± 0.076 (syst) for 
Am 2 3 = 2.4 x 10~ 3 eV 2 . 

The analysis is improved by comparing the positron 
spectrum in 18 variably sized energy bins between 
0.7 and 12.2 MeV to the expected number of D e events, 
again using Ato 2 3 = 2.4 x 10 -3 eV 2 . The analysis, per- 
formed with a standard x 2 estimator, uses four covari- 
ance matrices to include uncertainties in the antineu- 
trino signal, detector response, signal and background 
statistics, and background spectral shape. With very few 
positrons expected above 8 MeV, the region 8—12.2 MeV 
reduces the uncertainties in the correlated backgrounds 
with some additional contribution to the statistical un- 
certainty. 

The best fit results in 

sin 2 26»i 3 = 0.086 ± 0.041 (stat) ± 0.030 (syst) with 
a x 2 /DO F of 23.7/17, whereas the sin 2 26»i 3 = 0.0 
hypothesis results in a y 2 /DOF of 26.6/18. Using a 
frequentist approach [48j | we find an allowed region of 
0.017 < sin 2 26»i 3 < 0.16 at 90% CL, and exclude the no 
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Energy [MeV] 

FIG. 4. Top: Expected prompt energy spectra, including 
backgrounds, for the no-oscillation case and for the best fit 
sin 2 2#i3, superimposed on the measured spectrum. Inset: 
stacked histogram of backgrounds. Bottom: Difference be- 
tween data and the no-oscillation spectrum (data points) and 
difference between the best fit and no-oscillation expectations 
(curve) 



oscillation hypothesis at the 94.6% C.L. 

We determine our best estimate of the v e and back- 
ground rates with a pulls-based approach , the results 
of which are shown in Table Mil From the best fit we ob- 
tain a contribution from 9 Li reduced by ~19%, and with 
an uncertainty decreased from 52% to 26%. The fast 
neutron value is decreased by 5% with almost unchanged 
uncertainty. 



TABLE III. Summary of the effect of a pulls term approach 
on the fast neutron and 9 Li backgrounds and on the energy 
scale. Uncertainty values are in parentheses. 

Fast n. Bkg(%) a Li (%) EScale (value) 

Rate only tOO (46) 100 (52) 0.997 (0.007) 

Rate + Shape 95.2 (38) 81.5 (25.5) 0.998 (0.005) 



Figure [4] shows the measured positron spectrum super- 
imposed on the expected spectra for the no-oscillation 
hypothesis and for the best fit (including fitted back- 
grounds). 

Combining our result with the T2K |ll| and MI- 
NOS [IH measurements leads to 0.003 < sin 2 26»i 3 < 0.219 
at the 3<t level. 

In summary, Double Chooz has searched for 
v e disappearance using a 10 m 3 detector located 
1050 m from two reactors. A total of 4121 events 
were observed where 4344 ±165 were expected for no- 
oscillation, with a signal to background ratio of «11:1. 
In the context of neutrino oscillations, this deficit leads 
to sin 2 26»i 3 = 0.086 ± 0.041 (stat) ± 0.030 (syst), based 
on an analysis using rate and energy spectrum informa- 
tion. The no-oscillation hypothesis is ruled out at the 
94.6% C.L. Double Chooz continues to run, to reduce 
statistical and background systematic uncertainties. A 
near detector will soon lead to reduced reactor and de- 
tector systematic uncertainties and to an estimated la 
precision on sin 2 26>i3 of ~ 0.02. 
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